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Abstract. We present radio observations of the magnetic chemically peculiar Bp star 
HR Lup (HD 133880) at 647 and 277 MHz with the GMRT. At both frequencies the 
source is not detected but we are able to determine upper limits to the emission. The 
647 MHz limits are particularly useful, with a 5cr value of 0.45 mJy. Also, no large 
enhancements of the emission were seen. The non-detections, along with previously 
published higher frequency detections, provide evidence that an optically thick gyro- 
synchrotron model is the correct mechanism for the radio emission of HR Lup. 
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1. Introduction 

Sources of stellar radio emission can be split by the presence of non-thermal emission. Thermal 
emission is seen in massive s tars (Wolf-Ray et, O- and early B-type stars) that have strong ra- 



diatively driven stellar winds (I Wright 1 19751) . These winds are ionized and radiate via free-free 



emission. A number of early type stars also show a non-thermal component and this is usually 
due t o shock-acceleration associated with colliding stellar winds in binary systems (De Be cker! 



2007h . Several types of lower mass stars also show det ectable radio em ission, often non-thermal 



in na ture. Examples include the Sun, RS CVn st ars ('Slee et aL 2008 11). dMe stars ( Osten et al 



20061) . dwarf stars and brown dwarfs (IBerger et al..2006j . This emission is usually associated 
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with coronal emission and flare activity, is often time variable and also requires the presence of 
magnetic fields. 



Late B-type stars and early A-type stars have neither strong stellar winds nor substantial outer 
convective zones and so were not expected to be strong radio emitters. Magnetic chemically 
peculiar (MCP) stars have been known as radio emitters in the centimetre range since the mid 
1980s. MCP stars are a class of peculiar A and B stars (referred to as Ap or Bp stars) with 
strong (kGauss) magnetic fields. The class is characterized by large abundance anomalies in a 
range of elements which suggested that the photospheres of these stars are highly stable (with the 
turbulent motions stabilized by the magnetic field) allowing diffusion to take place. The origin of 
the magnetic fields in these stars is still unclear, with possibilities including either a fossil-field 
from the original collapse to form the star, o r possible dynamo action either in the radiative zone 
or at an earlier stage of evolution ('Arltl l2008l) . Although sub-surface convection zones in massive 
stars, associated with the iron opacity bump, could be associated with magnetic field gene ration, 
this mechanism is thought to not work for the late B-type stars considered here 



There is one MCP star that is of particular interest at radio wavelengths - CU Vir (AOp). This 
nea rby star has a perio d of 0.52 days and a strong magnetic field. From optical variabiUty tim- 
ing, IPvpereTaL ( 1998 ^ reported on a peculiar change in the optical (and presumably rotational) 
period. CU Vir also shows periodic, polarised outbursts as well as quiescent emission. The 
burst emission has been attributed to el ectron cyclotro n maser emission. Th e outbursts are see n 
over a wide range of frequencies (see (* Lo et al. 2012 ): ( Trigiho et al. 201 1 ): ( Ravi et al. 2010l) : 
(iStevens & George 2010) and references therein). 



There have been a small number of radio observations of late B and early A-stars. Drake et al. 



( 119870 observed 34 sources detecting only 5 at 5 GHz with the Very Large Array (VLA). They 
interpreted the observed radio emission as being due to gyro synchrotron emission from the mag- 
netosphere of the star due to a conti nuously inj ected population of mildly relativistic particles that 
are trapped in the magnetosphere. IWilson et a l. (1988) also using the VLA at 5 GHz surveyed 
16 M dw arf stars detecting o ne source, Gliese 735. A more recent VLA survey of MCP stars at 
5 GHz bv lLeone et al.l (119941) extended this to cover 40 stars, of which only 8 were detected. 



MCP stars on average show moderate circular polarisation, with radio properties similar to 
that of active cool stars such as RS CVns. They sh ow periodic changes in the detected magnetic 
field over the stellar rotation period (IBorra & Landstreet 1980) . The magnetic field topology of 
this kind of sta r is generally taken as a magnetic dipole tilted with respect to the rotation axis 
jBabcockl 19 49V though in the case of HR Lup the field is quadrupolar in nature (see below). In a 
few cases, MCP stars possess a n anisotropic stella r wind as a consequence of the wind interaction 
with the dipolar magnetic field dShore et al ]ll990l) . In about 25% of MCP stars non-thermal radio 
emission is observed and the rate of det ection seems to be correlated with the effective stellar 
temperature ( Linskv et al.lll992 : Leone et al.iil994) . 



The radio emission is generally interpreted as gyrosynchrotron emission from mildly relativ- 
istic electrons. These electrons are accelerated in current sheets, formed when the gas flow brakes 



Low -frequency GMRT observations of the magnetic Bp star HR Lap 



3 



the magnetic field lines, close to the magnetic equator and the electrons propagate along the mag- 
netic field lines towards the inner magnetospheric regions (Havnes & Goertz 1984). 



X-ray studie s have found that MCP stars are weak X-ray emitters with a detection rate of just 
over 10% (Drak e et al.lll994l) . The X-ray emission of MCP stars does not correlate with other 
stellar properties. Their radio properties also do not foll ow the Guedel-Benz s caling law as this 



would imply X-ray luminosities as high as 10^^ erg s ' ( Guedel & Bend 1993). The only X-ray 



detections of MCP stars are not particularly bright (L^ > lO^"-^ erg s lOrake et al]|2006h . 



In this paper we focus on low frequency observations of one specific object (HR Lup) which 
may give some insights into the processes going on in the star. 



2. HR Lup 



HR Lup (HD 133880, HR 5624; RA: 15 08 12.124, Dec: -40 35 02.15) is a rapidly rotating 
B-type chemically peculiar star of spectral type B8 Ivp. 

The star is rapidly rotating (u sin / ^ 103 kms"') and has a rotational period of 0.87 77 days, 
which is observed in optic al photometry, magnetic field measurements and radio flux (.Bailev et al. 
2012l;ISchmittet allbOOSb . 



The inferred steUar parameters for HR Lup are = 3.2M0, = l.OlRo, T,ff = 13000 K 
and age of 15.8 Myr (d erived from the fact that it is a member of the Upper Cen Lup association; 
Landstreet et al 1l2007h . HR Lup has a very strong magnetic field, typically 2.4 kG (iSchmitt et al 



20051) . This magnetic field varies from about 4 to +2 kG dLandstreet et al. Il990l) . UnUke most 
MCP stars the magnetic field is not dipolar, but quadrupolar (and this has c onsequences for th e 
observed radio variability). The rotational and magnetic axes are misaligned (iBailev et al. l2012h . 



HR Lup is a photometric variable with variations on the order of 0.15 mag in the U - 
band, which is p robably the result of the large magnetic field and surface abundance anomalies 
(IWaelkenslll985l) . 



HR Lup is a known r adio source , previ ously observed with the Australia Telescope Compact 



Array (ATCA). At 5 GHz lLim et al. I (119961) demonstrated that both the total intensity and circular 



polarisation of the source varied significantly and coherently according to the known rotational 
period. The total intensity varied from ~ 1 mJy to ~ 5 mJy with the degree of circular polarisa- 
tion reaching up to 20%. They noted that the emission shows broad peaks (suggesting a dipole 
contribution to the field) and narrower peaks at the predicted phases of a quadrupole contribution 
to the magnetic field. At 8 GHz the source is see n to have a flux d e nsity of 4.08 + 0.16 mJy 
though no indication of any variabil ity is presented ([Drake et al ]l2006h . iBailev et al. I (l2012h have 



re-reduced the ATCA radio data of iLim et al. 
the reduced data. 



(Il996h but there are no significant differences in 
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In this paper we present observations at 647 and 277 MHz with the intention of finding if 
there is any lower frequency emission and help to find if there is a cut-off" frequency. 



3. Observations and data reduction 

The Giant Metrewave Radio Telescope (GMRT) is located near Pune, India and consists of thirty 
45 m diameter radio dishes. The GMRT has a maximum and minimum baseline of 25 km and 
100 m, and has operating receivers at 150, 235, 325, 610 and 1400 MHz. 

We observed HR Lup with the GMRT on 2009 December 5th and 7th simultaneously at both 
647 and 277 MHz, implying only total intensity maps were constructed. At 647 MHz a bandwith 
of 16 MHz was used and at 277 MHz a bandwith of 6 MHz was used (both observations used 128 
channels). An integration time of 16 seconds was used. The total time on source was 4.67 hours 
over the two nights. During the observations 27 antennas were used. The observations consisted 
of flux density and bandpass calibration using 3C 286 at the start and end of the observations, and 
phase calibration using VSOP J1501-3918 every 30 minutes. 

Each spectral window was calibrated and imaged separately using the Common Astronomy 
Software Applications package (CASaQ). Any interference in the data was removed by manual 
inspection. Several rounds of phase only self-calibration were completed on the target data. 
During these iterations any visibility measurements that showed unusual phase excursions were 
rejected. The final image used for analysis below was primary beam corrected and has a central 
rms of 90 yuJy at 647 MHz and 6 mJy at 277 MHz. 

The source positions for the brighter sour ces were matched wi th that of the Sydney University 
Molonglo Sky Survey (SUMSS) at 843 MHz jMauch et al. Il2003h . No significant source position 



offset was found between the GMRT data and this survey. 

4. Results and Discussion 



At 647 MHz we do not detect the source. At the 5cr level this gives us an upper limit of the flux 
density of 0.45 mJy. Also, at 277 MHz we do not detect the source. At the 5cr level this gives 
us an upper limit of the flux of 30 mJy. The data were also cut into shorter time-scales (at 120 
seconds and the scan length ~ 600s) and imaged. No significant time variation was seen at the 
position of the source thus ruling out any bright bursts. 



Bailey et al. I (120 12h have produced the most recent and accurate ephemeris for HR Lup, 



combining several data sets to yield: 



* http://casa.nrao.edu/ 
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JDmin = 2445472.000(10) + 0.877476(9) ■ E, (1) 

where zero phase corresponds to minimum photometric brightness and also the minimum of the 
longitudinal magnetic field. With this ephemeris the principal radio maxima at GHz frequencies 
occur at phase = 0.0 and 0.5. 



From this, we can determine the phase of our observations. The observations are on two sep- 
arate days and they cover phases of 0.93-1.06 and 0.15-0.35 respectively. Thus the observations 
cover at lea st one of the period s of radio maxima previously observed. Using the 3cr period errors 
quoted by Bailey et al. I ( 120121) we estimate an uncertainty in the phase for these observations of 
around 0. 1 . 



Since we compared our positions to that of SUMSS it is worth noting that no source is given 
in their catalogue at the position of HR Lup, though investigation of the SUMSS mosaics indicates 
that there is a possible enhancement with a peak of 4.5 mJy/beam at the optical position of HR 
Lup. This source is below the detection threshold of 10 mJy/beam. 

If indeed there is a source corresponding to 4.5 mJy/beam at 843 MHz then this would imply 
an extremely steep cut-off between the 843 and 647 MHz observations (i.e. a > 7, for Syocv", 
which seems rather extreme). Though with no formal detection in the SUMMS maps this does 
not discount the possibility that the cut-off frequency is not located somewhere between 647 MHz 
and the previous detections at 5 GHz. If we discount the SUMSS point, then the implied limit on 
the spectral index between 647 MHz and the previously reported 5 GHz points is o- > 1 . These 
sources are possibly variable in time (although the ATCA observations show only a factor 2-3 
variability across the rotational period) over longer time-scales. 

The non-detection, particularly at 647 MHz, implies that the radio emission region intersect- 
ing the line of sight has a finite extent which is consistent with the optically thick gyrosynchrotron 
model of luinskv et al. (1992.) . The radio spectrum from this model will have a frequency (vpeak) 
where the emission peaks, the location of which depends on the wind density/geometry and the 
magnetic field properties. The radio spectrum then falls away on either side of this peak fre- 
quency, and the spectral slope on either side of the Vp^ak with the electron energy spectrum. It 
is likely that Vp^ak lies somewhere between 647 MHz and 5 GHz, and observations with a much 
broader bandwidth will be important to fully constrain the detailed shape of the emission. 



5. Summary 



In summary, we present non-detections of the MCP star, HR Lup with the GMRT at 647 MHz 
and 277 MHz. No significant time variation was seen at the position of the source thus ruling out 
bright bursts, as seen in CU Vir, in this dataset. We consider that these non-detections provide 
evidence that the emission mechanism for this star is optically thick gyrosynchrotron emission. 
We suspect that the peak emission frequency lies between 647 MHz and 5 GHz. Follow-up ob- 
servations at frequencies in this range would be particularly useful to constrain the spectral model 
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Table 1. The GMRT observations of HR Lup in UT at 647/277 MHz. Given are the start and end times of 
the observing and the range of phase covered by these observations, using the ephemeris of iBailev et al. I 
l20l2). 



Start Observing 


End Observing 


Phase range 


04:35:10.9 05 Dec 2009 
03:23:28.3 07 Dec 2009 


07:22:40.4 05 Dec 2009 
07:35:41.3 07 Dec 2009 


0.93-1.06 
0.15-0.35 



and indeed provide constraints on the energy injection of electrons in the stellar magnetosphere. 
Broad-band observations (covering both the GHz and sub-GHz regimes) are necessary to con- 
strain the emission mechanisms in MCP stars and coverage of the entire rotational period are 
necessary to detect the presence of short-Uved intense bursts, which may well be a feature of the 
radio emission from the strongly magnetic stars. 
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